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ABSTRACT 

The shock model of gamma-ray bursts (GRBs) contains two equipartition parameters: the magnetic 
energy density and the kinetic energy density of the electrons relative to the total energy density of 
the shock, es and te, respectively. These are free parameters within the model. Whereas the Weibel 
shock theory and numerical simulations fix es at the level of ~fewx (10~'^...10~'*), no understanding of 
Ef, exists so far. Here we demonstrate that it inevitably follows from the theory that Ce ~ The 
GRB afteglow data fully agree with this theoretical prediction. Our result explains why the electrons 
are close to equipartition in GRBs. The eg — relation can potentially be used to reduce the number 
of free parameters in afterglow models. 

Subject headings: gamma rays: bursts — shock waves — magnetic fields 



1. INTRODUCTION 

Shocks in gamma-ray bursters (GRBs) as well as in 
many other astrophysical sources are collisionless. Their 
physics is much more complex than that of standard hy- 
dro dynamic shocks. T he Weibel shock theory developed 
by iMedvedev fc Loehl (^999) has been confirmed by a 
large number of nun ierical simulations both in the ultra- 
relatiyistic regime llSilya. et al.l l2003t iNishikawa. et alJ 
l2003t iFrederiksen. et all '2004^ 'Spitkovskv '2005') and in 
the n on-relativistic one (Medvedev. Silva. Kamionkowski: 
l2006j) . In particular, it has been predicted and then nu- 
merically confirmed that the magnetic equipartition pa- 
rameter, es, defined as the ratio of the magnetic energy 
density to the total kinetic energy density of a shock, 



TTipC^nT ' 



(1) 



falls in the range of '--^fewx (10"'^... 10"*), where n is the co- 
moving particle density and T is the shock Lorentz factor; 
we assumed F 3> 1. 

Multiwavelength spectral fits to afterglows allowed 
reliable determination o f micro-physics parameters for 
a number of GRB s ("Panaitescu & Kumar" '2001a"b', 
2002: Chc v alicr"FLil ITOO O: Li feChevalier 2001, 2003^ 
Yost, et ahl l2003t iFrail et alJ l200i iBereerT et all i2004i 
McMahon. et all 12004 lPanaitescull2005j) . The analvses 
found that cb is indeed falls in the range ^ 10^^. ..10^^, 
however with relatively large errorbars — typically of 
about a decade, depending on quality of data (afterglow 
temporal and spectral coverage) and the method of analy- 
sis. These results confirm (within uncertainty) the predic- 
tions of the Weibel shock theory. We note that this is the 
only theory, which explains the origin of magnetic fields in 
relativistic shocks. 



The afterglow analyses also yield the values of Cg, de- 
fined as the ratio of the kinetic energy density of electrons 
to the total energy density of the shock: ^ 



rripC^nT 



(2) 



The values turn 
-fewx (10-1. ..10^2)^ 



out to be clustered at 
Given the fact that before being 



shocked the electrons carry only about ^ rrie/mp ~ 10" 
of the total energy of the ejecta, they must be acceler- 
ated somehow in the downstream region. Traditionally, 
one invokes diffusive Fermi acceleration for this purpose. 
However, the Fermi theory cannot accurately predict the 
value of Cp in GRB shock s. Moreover, recent studies 
ijNiemiec fc Ostrowskill2006t iLi fc Waxm'anll2006() put the 
efficiency, or even the very presence, of Fermi accelera- 
tion at GRB shocks into question. Thus, until now, we 
had little understanding of the electron acceleration in 
GRBs and we did not know why the electrons are close to 
equipartition. 

In this paper, we demonstrate that eg — -^v^es (with the 
constant A ~ 1) in relativistic, baryon-dominated shocks 
(i.e., shocks in electron-proton plasma). The result is very 
robust and is based solely on well-known properties of col- 
lisionless shocks. 

2. EVALUATION OF 

Magnetic fields are generated at shocks by the Weibel 
instability. In baryon-dominated shocks, the value ~ 
10"'^ is limited by the charge-separation effects, which 
mo d ify the instability growth ra t e and the its dynam - 
ics (jWiersma fc Achtenber j l2004t iTzonfras. et al.ll2006j) . 
Thus, saturation occurs at the equipartition with the light- 
est species — the electrons. At such low fields, protons 
keep streaming in current filaments,'^ whereas the elec- 



^Also: Institute for Nuclear Fusion, RRC "Kurchatov Institute", Moscow 123182, Russia 

^Often, te is defined for a power-law distributed electrons, so the power-law index enters the definition: 7e,min = 
((p - 2)/(p - 1)) (mp/me)ree, where 

^e^min IS the low-energy cut-off of the electron distribution and p is the power-law index. In this case, 
to measure the total energetics one uses another (related) parameter e^. Here we make no assumption about the electron energy distribution, 
hence our ee measures the overall energetics of the electrons. Thus our te, defined in Eq. 0, is identical to used by some authors. 

^By the way, this explains why numerical simulations still cannot ful ly r esolve baryon-dominated shocks, in contrast to the electron-positron 
pair-dominated shocks (see, e.g.. iFrederiksen. et al...2004 : .Spitkovskv.i.2005l) . 
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trons, being much lighter than the protons, are quickly 
isotropized in the random fields and form a uniform back- 
ground. The current filaments are formed by the protons 
moving roughly at the speed of light (their Lorentz factor 
is r). Hence, the y are sources of both the magnetic and 
electr ostatic fields (|Hededal. et aUlToO^ iNishikawa. et al.1 
l20nl . These fields are related to each other as 

B = PE, (3) 

where fi — \/l — T^^ ^ 1. An electron, moving toward a 
filament gains energy 

Ue ~ elE ~ elB. (4) 

The typical radial distance the electron travels is about 
half the distance between the filaments, I ~ X{c/ujpp^rei)i 
where A ^ 1 is the dimensionless parameter, c/ujpp^rei 
is the relativistic proton skin depth — the typical scale 
of structures in the Weibel turbulence, and ojpp.rei — 
(47re^n/mp7p)^/^ is the relativistic proton plasma fre- 
quency and 7p ~ r. The parameter A accounts for the 
actual geometry of the filaments, the electrostatic shield- 
ing in plasmas, the effects of the electrons on the current 
distribution, etc. All these effects introduce only a small, 
factor of two, uncertainty, as is discussed at the end of 
this section. Finally, the electron energy density behind 
the shock front is 

Ue — nue — XeBnc/upp^rei- (5) 
This equation can be cast into the form 

Ce ^ Xy/eB- (6) 

Note that we didn't make any assumptions here on 
whether the shock compression have already occurred or 
not (i.e., how far downstream we are). We just used the 
fact that the shock magnetic fields are due to proton cur- 
rents, which also produce electrostatic fields. These elec- 
trostatic fields locally accelerate electrons on their way 
in and decelerate as they go away from the filament. 
Since the electron emissivity is F^, (x B^j^, the electrons 
strongly radiate near the filaments, where their energy and 
the magnetic field are both at maximum. Hence, Eq. (jSJ 
represents the emission-weighted relation between and 
es . It is this relation that should be found in GRB obser- 
vations. 

In the electron-positron-dominated shocks, the situa- 
tion is drastically different. By symmetry (e"*" and e~ 
have identical masses), the current filaments are formed 
by nearly equal numbers of electrons and positrons mov- 
ing in opposite directions. Hence, the net charge of such 
filaments is vanishing and no substantial particle accelera- 
tion is expected. T his is in full agreement with simulations 
llSDitkovskvll2n05D . 

A rather important question is: How well is the pa- 
rameted A constrained? We performed a special study 
using a more detailed model of electron dynamics in cur- 
rent filament fields. The details of this work will be re- 
ported elsewhere. Here we just cite the relevant result: 
the emission-weighted value of A is very insensitive to the 
details, e.g., the current distribution, filament filling fac- 
tor, etc. For typical parameters obtained from PIC simu- 
lations, A ranges between 1 ^ A ^ 3 while the filling factor 



ranges from ~ 1/4 (near the shock front) to ^ 1/100 (far 
downstream). Thus, A is constrained rather well, it is not 
just a new free parameter of a model. 

3. COMPARISON WITH OBSERVATIONS 

To compare with observations, we have taken the most 
recent and the best analyzed sample of data containing 
ten afterglows (Panaitescu 2005). These GRB afterglows 
were fitted to a number of afterglow models, which include 
combinations of the models of the external media profiles 
(constant density and wind-like) with the models for the 
ejecta structure (jetted, structured outflows and energy in- 
jection models). For each model, the reduced was given 
and for models with "reasonably good fits" (x^/dof< 4), 
the micro-physical parameters are given. 

From this data set, we have chosen the best fit model 
(having the smallest reduced value) for each GRB. The 
parameters are given in Table Note that for two GRBs, 
there are two equally good fits, hence we included both. 
As one can see, es varies over two orders of magnitude and 
ee (which is in definitions of the paper Panaitescu 2005} 
varies by one order of magnitude. If es and eg are sta- 
tistically independent, the scatter of the quantity \f^B 
should be of about two decades. We plot e^l for the 
best fit models in Figure QJi. Clearly, the data is clustered 
around unity with only little scatter. No indvidual confi- 
de nce intervals for t he best fit parameters were provided 
in iPanaitescul l)2005() . Only the overall uncertainties were 
provided: ailges) = 1, cr(ee) — 0.3ee. These yield an 
uncertainly of the ratio a{\g[ee/ ^/Tb]) — 0.5. 

The clustering of ee/y^eB near unity is neither acciden- 
tal, not an artifact of fitting. To demonstrate this, we plot 
in Figure [TJd the values of tel \f^B for all models reported 
by Panaitescu (200l|). The goodness of the fit for these 
models is x^/dof< 4. The data points are scattered over 
almost three decades, which is consistent with degradation 
of statistical correlation of tf. and in poor fits. 

Finally, we performed a linear fit in the log eg — log e b- 
space. We set the intercept to zero (log A ~ 0) to re- 
duce the number of degrees of freedom, because a two- 
parameter fit does not give a statistically acceptible result. 
The one-parameter fit yields the exponent in the relation 
ee = t% being s = 0.49 ± 0.07 with the p-value - 10"'^. 
This is in excellent agreement with Eq. @. 

4. CONCLUSIONS 

In this paper, we have shown that in relativistic, baryon- 
dominated shocks te — X^/Fb with the constant A ~ 1. 
The result inevitably follows from the micro-structure of 
collisionless shocks. No ad hoc assumptions were made. 
Interestingly, the values of t^j \/^b derived from afterglow 
data of ten GRBs are clustered around unity, thus sup- 
porting our theory. Since the typical value of es is about 
10"^, the corresponding value of eg should be ~ 0.03. In- 
deed, the afterglow fits show that 0.03 is the typical value 
of the Ce parameter. Thus, our theory explains why the 
electrons in GRB shocks are close to equipartition. In- 
terestingly, the typical uncertainly in determined from 
afterglow fits is often rather large, - of about an order of 
magnitude or even more. In contrast, the theoretical un- 
certainly in the value of A is only a factor of few at most, 



3 



and is very likely even less. Therefore, one can use the 
obtained relation ^ -^/es to reduce the number of free 
parameters in afterglow models. It will be very interest- 
ing to investigate how the goodness of afterglow model fits 
changes when the above relation is used. 

A number of important questions are left outside of the 
present study. It would be interesting to calculate the en- 
ergy distribution of electrons. Some simple estimates can, 
however, be given as follows. Far away from the filaments, 
the electrons form the isotropic distribution with a typ- 
ical Lorentz factor je,mm ~ T, because they have been 
pitch-angle scattered in small-scale magnetic fields (these 
scatterings do not change particle energy). Inside the fil- 
aments the 7-factor of the electrons shall be maximum, 
^e,max ~ {mp/me)T^£e- Since these electrons contribute 
the most to the observed emission, the emission-weighted 
(that is, "observed") distribution will be peaked at 'ye,max, 
which shall correspond to Epeak in the spectral distribu- 
tion of GRB emission. Other mechanisms are likely needed 
in order to explain the po wer-law spectra abo ve Epeak ■ Nu- 
merical PIC simulations l|Silva. et alJl2003l indicate that 
reconnection events occuring during current filament coa- 
lescences lead to local acceleration of electrons. This can 



possibly produce a power-law distribution of electrons. Al- 
ternatively, a power-law radiation spectrum can be pro- 
duced even by a monoenergetic electron distribution, pro- 
vided the spati al spectrum of sm all-scale magnetic fields 
is a power-law llMedvedevI l2006j) . PIC simulations indi- 
cate that power-law s pectra of magnetic felds indeed form 
at relativistic shocks (jFrederiksen. et al]l2004j) . A related 
question is how large the fraction of the energetic elec- 
trons is. A simple estimate indicates that is should be 
of order the filling factor of current filaments, which can 
be determined from PIC simulations: it varies with the 
downstream distance from ^ 1/ few to ^ 1/100 in typical 
electron-proton runs. However, more accurate answers to 
the above questions require detailed modeling of the elec- 
tron dynamics, taking into account the residence time of 
electrons inside the filaments (the electrons are defiected 
in the direction of the ion current as they move toward 
a filament). We are developing a detailed model and the 
results will be presented in forthcoming publications. 

This work has been supported by NASA grant NNG- 
04GM41G and DoE grant DE-FG02-04ER54790. 
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Table 1 

I2005D . The models selected by the lowest reduced ^ (fourth column). 



The best fit afterglow model parameters from lPanaitescu| 
Two GRBs, 991216 and 000926, have more than one best fit model. The models in the last column are SO=structured outflow, 
EI=model with energy injection, JET=model with constant F within a jet opening angle, ISM=constant density medium, 

WIND=wind-like medium. 
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Fig. 1. — The ta/y/tB ratio for ten GRB afterglows analysed bv IPanaitescul 120nfil1 . (a) — The parameters of the best fit model (with 
the lowest x'^/dof) are used. The clustering of data points around unity is evident. All points are consistent with unity within the quoted 
errorbars. (b) — All reported models are used (with X'^/<iof< 4). The scatter is substantially larger, indicating the degradation of the Ce — es 
correlation for poor fits. 



